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Abstract 



We review recent measurements of the high-frequency dynamic magnetic 
susceptibility in the high-T c superconducting systems La2- :r Sr, E Cu04 and 
YBa2Cu306+iE- Experiments were performed using the chopper spectrom- 
eters HET and MARI at the ISIS spallation source. We have placed our 
measurements on an absolute intensity scale, this allows systematic trends to 
be seen and comparisons with theory to be made. We find that the insulating 
S = ^ antiferromagnetic parent compounds show a dramatic renormalization 
of the spin wave intensity. The effect of doping on the response is to cause 
broadenings in wave vector and large redistributions of spectral weight in the 
frequency spectrum. 

I. INTRODUCTION 

The collective excitations of the spins in the CuO"2 planes of the high temperature su- 
perconductors display a plethora of behavior [I], perhaps as exotic as the bulk properties 
themselves. A characterization of these excitations is motivated by several considerations. 
The spin excitations, which can be measured by a well understood and easily interpretable 
probe such as neutron scattering, provide a window on the electronic correlations in these 
complex materials. Indeed, Cooper pairing itself is a correlation between electron spins. A 
further motivation, especially in view of the unconventional pairing in these materials, is 
that the magnetic excitations may be involved in the pairing attraction. 

In the present paper we review our recent measurements of the dynamic response 
over a wide frequency range i.e. up to energies of order 2 J. We find that in the supercon- 
ductors there is significant spectral weight above the superconducting gap energy. For both 
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superconductors studied, the magnetic response x"(q,w) is actually stronger for frequencies 
just above 2 A than in the parent insulating antiferromagnet. We have been careful to con- 
vert our measurements into absolute units. This allows a systematic comparison between 
different systems, compositions and theories to be made. 



II. EXPERIMENTAL 

Our experiments were performed using the HET and MARI spectrometers at the ISIS 
pulsed spallation neutron source of the Rutherford Appleton Laboratory. HET and MARI 
are direct geometry chopper spectrometers. A pulse of neutrons of approximately 1 [is dura- 
tion is produced when a 800 MeV pulsed proton beam hits a tantalum target. The neutrons 
are monochromated using a Fermi chopper at 10 m from the target and appropriately phased 
to the proton pulse. Scattered neutrons are detected using 800 3 He detectors 4 m from the 
sample. 

Neutron scattering directly measures the imaginary part of the generalized magnetic 
susceptibility x"(Q,u;). The scattering cross-section for an isotropic system is (7], 



where (7r e ) 2 =0.2905 barn sr _1 /i^ 2 , and k/ are the incident and final neutron wavevec- 
tors and |F(Q)| 2 is the magnetic form factor. Absolute unit conventions were performed by 
measuring the elastic incoherent scattering from a vanadium standard || under the same 
experimental conditions and by measuring the low-frequency coherent phonon scattering 
from the sample ||. Throughout this paper, we label momentum transfers (Q x , Q y ,Q z ) 
in units of A -1 by their reciprocal space positions Q=(h,k,l)= (2nQ x /a,2TcQ y /b,TcQ z /c). 
Following previous practice, we use the orthorhombic nomenclature to label reciprocal 
space in the La2- x Sr x Cu04 system and tetragonal nomenclature in the YBa 2 Cu306+a; sys- 
tem. This means that the Cu0 2 planes are parallel to (010) in La 2 _ x Sr x Cu04 and (001) 
in YBa 2 Cu 3 6+:1 .. Details of the samples are given elsewhere [|3|-fJ. In the case of the 
La 2 _ a; Sr 3 .Cu04 system, data were collected with the (001) plane coincident with the princi- 
pal scattering plane of the spectrometer. For YBa 2 Cus06+x, the (110) plane was used. 



III. THE QUANTUM ANTIFERROMAGNETS 

We first discuss the parent antiferromagnets La 2 Cu0 4 and YBa 2 Cu 3 6 .i5. F° r the pur- 
pose of this paper, we model the systems as a set of weakly coupled Cu0 2 layers or bilayers. 
In this case, the high-frequency spin excitations can be described by the Heisenberg Hamil- 
tonian for a single Cu0 2 layer or bilayer, 

H = J £j\\S i -S j + J £j ± S i -S jl . (2) 

ij ij' 

The first term in Eq. |2| represents the nearest-neighbor coupling between Cu spins in the 
same Cu0 2 plane. The second term (not present in the single layer compound) represents 
the coupling between nearest-neighbor Cu spins in different layers. 



2 



In the case of the single layer compound (J± = 0), conventional spin-wave theory in the 
classical large-5 limit yields the transverse dynamic susceptibility 

1 f'2 

X '[(Q,uj) = Z x ^g^lS ( 1^^ ) 8(fiu>±hw{Q)), (3) 



where, 

r i 1/2 

hu(Q) = 2Z C J\\ 1-7 2 (Q) , (4) 



and 7(Q) = cos(7r/i) cos(7r/). We have included a "quantum renormalization" of the overall 
amplitude Z x . In the conventional linear spin- wave model applicable for large 5, Z x = 1. 
For small 5, quantum corrections [[T0| , [TT]1 become important (the Neel state is not a good 



approximation to the ground state) leading to a renormalization of the overall scale of the 
spin-wave dispersion and to a reduction in magnetic response with respect to the classical 
spin-wave theory. The renormalization of the overall scale can be included in the exchange 
constant J* = Z C J, where J is the exchange constant occurring in Eq. |^. The value of Z x 
can be obtained by neutron scattering, if measurements are placed on an absolute intensity 
scale. However, Z c cannot be measured directly from inelastic neutron scattering and must 
be estimated from theory. In the case of the 5 = | square- lattice antiferro magnet, Singh 
|T0| and Igarashi |L1] have estimated Z c = 1.18 and Z x =0.51 based on a 1/5 expansion. 

The presence of the second term (i.e. J± ^ 0) in Eq. |] leads to the existence of two 
branches in the spin-wave dispersion which can be labeled according to whether neighboring 
spins in different planes rotate in the same direction ("acoustic" or "odd" mode) or in op- 
posite directions ( "optical" or "even" mode) about their time-averaged (ordered) directions. 
In the conventional linear spin-wave approximation, the acoustic and optic modes have the 
response functions per formula unit |I2[ 



Xac(Q, o) = Z x 7T g fi B S I + 7 (q) ) ( 5 ) 

( -k Azl\ 

xsin 2 ] 5(hu±huj ac {Q)) (6) 



and 



fn Azl\ 

x cos 2 - j 5(hu±hcu op (Q)) , (8) 
respectively. The dispersion relations are 

^ac(Q) = 2Z C J|, {l - 7 2 (Q) + J ± /J,| [1 ± 7 (Q)]} 1/2 , (9) 

where 7(Q) = |[cos(27r/i) + cos(27rfc)] and Az= 3.2 A is the separation of the CuC>2 planes 
in a bilayer. The inter-planar coupling term in Eq. ^| leads to no additional dispersion along 
the z-direction, only a modulation in the amplitude of the response which nevertheless can 
be used to distinguish between the two modes. 
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A. La 2 Cu0 4 



The two-dimensionality of the scattering in La2Cu04 means that we are able to cut 
through the spin waves at several energy transfers for a single spectrometer setting. Fig. [I] 
shows data collected for E,i =300 meV and kj || (010) on the MARI spectrometer. Panels (b)- 
(f) show constant energy cuts along the (1,0,0) direction. A spin-wave peak is observed near 
h = 1. Twin peaks due to spin waves propagating in opposite directions are not observed 
due to the poor out-of-plane resolution in the (001) direction. The broadening in the peak at 
higher frequencies is due to the spin-wave dispersion. A convenient way to display the data in 
Fig. [I] is as a local- or wavevector-integrated susceptibility = J BZ x"(Q> u ) d 3 Q/ J d 3 Q, 

this is shown in Fig. 0(a). The solid lines in Fig. |l](b)-(f) and Fig. ^| are fits to the spin-wave 
model described by Eq. |3|. Quantitative analysis yields an exchange constant consistent with 
our previous determination || J* = 156 ± 5 meV and Z x = 0.39 ±0.1. This is in agreement 
with the calculation of Igarashi |TT| and demonstrates the importance of quantum corrections 
in this system. 



B. YBa 2 Cu 3 6 15 

As discussed above, the bilayer nature of this material results in two collective modes with 
additional structure to x"{Q,i u ) as compared with the single layer material. The optic mode 
has a gap at the 2-D magnetic zone centers such as of huj g = 2JJ* L J*\. Its intensity dis- 
plays an overall modulation cos 2 (7r Azl/c) which has maxima at positions I = 0, 3.7, 7.3, . . .. 
Correspondingly, the acoustic mode has an overall modulation of sin 2 (n Azl/c) and there- 
fore has maxima at I = 1.8, 5.5,.... We are able to probe the spin waves at various values 
of hu and I by varying the incident energy. Fig. ^|(b)-(e) show data collected for different 
energy transfers at I values corresponding to positions where the scattering is predominately 
acoustic. As in La2Cu04, a peak is observed near the magnetic zone center, for all energy 
transfers investigated, due to propagating spin waves. Similar cuts probing I values where 
the optic mode dominates, are shown in Fig. |||(g)-(j). In this case, no peak is observed at 
the lowest energies because these energies are below the optic gap. A more convenient way 
to look at these data is in the form of a local susceptibility (see above). Fig. |4] shows the 
contributions of the two branches to the local susceptibility extracted from the data in Fig. |3|. 
When the data are plotted in this way, we see that the optic and acoustic contributions are 
equal at higher frequencies. While below hu g , the contribution from the optical mode is 
zero within the experimental error. A detailed analysis yields a value for the optic gap 
of hujg = 74 ± 5 meV. A similar value for the optic gap was obtained using reactor-based 
instrumentation [ |i~3|| . 

In order to determine the exchange constant J]_ describing the coupling of spins in 
different planes, we also need to know Jjj 1 because the optic gap is Tiuj g = 2^JJj_J*\. The in- 
plane exchange coupling, J,f can be obtained from a measurement of the high-frequency spin 
waves. Fig. |5| shows data collected with the higher incident energy Ei =600 meV. We note 
that for energy transfers above 245 meV there is little variation in the intensity, suggesting 
that this energy is close to the zone boundary energy. A simultaneous resolution-corrected 
fit of the linear spin wave model (Eqs. @-(0)) to all our data yields values for the exchange 
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constants of Jm = 125 ± 5 meV and J]_ = 11 ± 2 meV. Further, since our measurements are 
in absolute units we are able to estimate the amplitude renormalization. We find a value 
of Z x = 0.4 ± 0.1, as for La2Cu04, again the quantum renormalization is close to the value 
predicted by a 1/S expansion [Z x = 0.51). 



IV. THE SUPERCONDUCTORS 



Of particular interest is the nature of the spin fluctuations for metallic and super- 

conducting compositions . We have therefore studied the compounds Lai. 86 Sr .i4CuO4 (a 
composition close to optimal doping) and YBa 2 Cu30 66 (an underdoped superconductor). 
In both cases we find strong high-frequency spin fluctuations in the superconducting state. 



A. Lai. 86 Sro.i 4 Cu0 4 



Fig. |l|(h)-(l) shows data collected for Lai.86Sro.i4CuO"4 under the same conditions as 
Fig. [TJ and in the same units. Doping has a dramatic effect on the magnetic excitations: at 
low frequencies, the peaks are broader in reciprocal space and at higher frequencies there 
is a large suppression of the intensity. At the highest frequencies, the peaks appear to 
be disappearing more rapidly in the superconductor than the insulator. This can be seen 
clearly when we extract and plot the local susceptibility in Fig. @(b). When we include low 
frequency data from a reactor-based experiment, it becomes clear that the local response is 
peaked near 22 meV. Thus doping leads to a shift of spectral weight to intermediate energies. 



B. YBa 2 Cu 3 6 . 6 



The low frequency dynamics of superconducting YBa 2 Cu 3 6+x have been the subject of 
considerable investigation [l|]. The YBa 2 Cu30 66 (T c = 62.7 K) sample used in the current 
study is of high quality: it shows a sharp "resonance peak" [HJ at 34 meV on entering the 
superconducting state and has recently been shown |16] to develop incommensurate peaks in 
for frequencies near 25 meV. We concentrate here on the nature of the dynamics for 
temperatures just above the superconducting transition temperature. Fig. ^| shows inelastic 
data collected for T = 80 K on the MARI spectrometer in a similar manner to that in Fig. |3|. 
The reader should note that the energies chosen for Fig. ^| are different from those chosen 
for Fig. ^ because of the different detector positions on the HET and MARI spectrometers. 
By assuming that the scattering has a Gaussian wavevector dependence for each energy, we 
have extracted the frequency variation of the acoustic and optical contributions to the local 
susceptibility x"( u ) from data such as those shown in the figure. Fig. |7| shows the results. 

Comparison of Figs. |6] and [7] with Figs. |3] and f| shows that the magnetic response has 
changed dramatically between the metal and the insulator. This is well-known at low fre- 
quencies |TI], as in La 2 _ x Sra;Cu04, there is a large broadening of the low frequency peak in 
wavevector. More dramatic is the redistribution of spectral weight in energy revealed by the 
local susceptibility in Fig. [j]. Inspection of Fig. |] and Fig. [7|, which are in the same units and 
therefore can be directly compared, shows that doping causes a shift of spectral weight from 
high frequencies (above 100 meV) into the range around 50-100 meV. The gap in the optic 
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fluctuations remains and is possibly slightly reduced. Reactor-based measurements yield a 
similar value for the optic gap [17 . 



V. DISCUSSION 

Our measurements of x"{ u ) 011 La 2 Cu04 and YBa 2 Cus06.i5 are summarized in Fig. ||(a) 
and Fig. f|. We note that observed x"( u ) f° r La 2 Cu0 4 and the acoustic branch of 
YBa 2 Cu30 6 .i5 are of approximately equal intensity: the units of the susceptibility are ex- 
pressed per formula unit in both cases, YBa 2 Cu 3 6+a; has two Cu in the Cu0 2 planes 
per formula unit. When we compare our measurements on the quantum antiferromagnets 
La 2 Cu0 4 and YBa 2 Cu 3 6 .i 5 , we find that, in both cases, linear spin-wave theory accounts 
well for the dispersion and amplitude variation. However, an overall renormalization of the 
amplitude Z x must be included to account for the observed intensity of the one-magnon 
response. Our observations agree with calculations of Z x =0.51 based on a 1/5* expansion. 
Although we are able to explain the one-magnon excitation spectrum, the 1/S expansion 



also predicts additional spectral weight ]TIJ which has not yet been observed. 

We now turn the metals Lai.86Sr .i4CuO4 and YBa 2 Cu306.6- The response of the copper 
spins obeys the sum rule, 



m 



3h r°° x!'(w)du 



7r J-00 1 — exp(— huo/kT) 



(10) 



where (m 2 ) is the mean squared moment on the Cu sites. In both superconductors studied 
here, doping and the consequential metal-insulator transition led to a large redistribution of 
spectral weight and formation of a peak in x"^) a ^ intermediate frequency. In the light of 
the sum rule, the existence of the peak in itself is not unexpected since spectral weight is lost 
from the Bragg peak present in the insulator and the high frequency response is strongly 
suppressed. What is surprising is the relatively small energy scale over which the response is 
distributed compared with other paramagnetic metals. Such small electronic energy scales 
are, of course, characteristic of the cuprate superconductors. 

Work at ORNL is supported by US-DOE under Contract No. DE-AC05-96OR22464 with 
Lockheed Martin Research, Inc. We are grateful for the finacial support of the UK-EPSRC 
and NATO. 
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FIGURES 



FIG. 1. (a)-(f) Magnetic scattering from La2Cu04 and (g)-(l) from Lai.86Sro.i4Cu04 Q. All 
scans are in the same absolute units. Note the magnetic peak is broader in the superconductor 
and suppressed at higher frequencies. 

FIG. 2. Local susceptibility § derived from data such as that in Fig. |l| (closed circles) and from 
reactor-based measurements Note huge redistribution of spectral weight caused by doping. 

FIG. 3. Constant energy scans showing magnetic scattering from YBa2Cu30g.i5 for wavevectors 
Q = (h,h,l). (b)-(e) I chosen to emphasize acoustic modes, (g)-(j) I chosen to emphasize optic 
modes. The onset of scattering at the optic position is about 74 meV. 

FIG. 4. Local or wavevector-integrated energy-dependent magnetic susceptibility at acoustic 
and optic /-positions in YBa2Cu306.i5- No integration over I has been performed. Points have been 
obtained by integrating over the spin-wave peaks and correcting for the Cu 2+ magnetic form factor, 
Bose factor and instrumental resolution. Solid lines are a fit to the spin wave model described in 
the text. 

FIG. 5. (a) Dispersion relation of YBa2Cu30g.i5. Closed circles and solid line are acoustic 
mode. Open circles and dashed line are optic mode, (b)-(e) Constant energy scans showing the 
high-frequency magnetic scattering from YBa2Cu30g.i5. 

FIG. 6. Constant energy scans showing magnetic scattering from YBa2Cu30g.6 for wavevectors 
Q = (h,h,l). (a)-(c) I chosen to emphasize acoustic modes, (d)-(f) I chosen to emphasize optic 
modes. 

FIG. 7. Acoustic and optic contributions to the local susceptibility (see Fig. |j) in YBa2Cu306.6- 
Solid lines are a guide to the eye. 
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TABLES 



TABLE I. The table summarizes experimentally-determined spin-wave parameters for 
La2- x Sr x Cu04 and YE^CusOe+x- In the case of the superconductors, J|* and Z c are determined 
from the highest frequencies studied. 

Jj_ TlUJg 

La 2 Cu0 4 156±5meV 0.39 ±0.1 

Lai. 86 Sr .i4CuO 4 130 ± 5 me V 0.15 ±0.06 

YBa 2 Cu 3 6 .i5 125±5meV 11 ± 2meV 0.4 ±0.1 74 ± 5 meV 

YBa 2 Cu 3 6 . 6 125 ±20meV 0.24 ±0.12 ~ 60 meV 
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